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The Goldberger-Treiman discrepancy in SU(3) is analyzed in the frame- 
work of heavy baryon chiral perturbation theory (HBChPT). It is shown that 
the discrepancy at leading order is entirely given by counterterms from the 



0(p ) Lagrangian, and that the first subleading corrections are suppressed by 



two powers in the HBChPT expansion. These subleading corrections include 

X' 

$_i ' meson- loop contributions as well as counterterms from the 0(p ) Lagrangian. 



Some one-loop contributions are calculated and found to be small. Using the 
three discrepancies (jrNN, KNA and KNU) which can be extracted from 
existing experimental data, we find that the HBChPT calculation favors the 
smaller g^NN values obtained in recent partial wave analyses. 
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I. INTRODUCTION 



The Goldberger-Treiman relation (GTR) obtained from matrix elements of the di- 
vergence of axial currents between spin 1/2 baryons, is an important indicator of explicit 
chiral symmetry breaking by the quark masses. It interrelates baryon masses, axial vector 
couplings, the baryon-pseudoscalar meson (Goldstone boson = GB) couplings and the GB 
decay constants. Explicit chiral symmetry breaking leads to a departure from the GTR 
(defined below) which is called the Goldberger-Treiman discrepancy (GTD). 

The GTD has been repeatedly discussed over time and for several reasons there 
were difficulties in arriving at a clear understanding. On one hand, there was no available 
effective theory with a systematic expansion to address the problem, and on the other hand 
the experimental values of the baryon-GB couplings were too poorly known. In recent years, 
progress has been made on both fronts. There is now a baryon chiral effective theory that 
permits a consistent expansion of the discrepancy 0-§f. There has also been progress in 
the determinations of the baryon-GB couplings that are the main source of uncertainty in 
the phenomenological extraction of the discrepancies. In fact, the current knowledge of the 
couplings g n NN, gKNA and, to a lesser extent, Qkny* is good enough to justify a new look 
at the GTD in SU(3). In this work we study the GTD in the light of heavy baryon chiral 
perturbation theory (HBChPT) J§§. 

Let us first briefly review the derivation of the GTR |7j and the definition of the GTD. 
We consider the matrix elements of the octet axial current = |g(x)7 A1 7s A a q(x) (the 
Gell-Mann matrices are normalized to Tr(A a A b ) = 2S ab ) between states of the baryon octet: 

(&, Pb | a; I a, Pa ) = u b ( Pb ) [~ i,gt(q 2 ) - ^gtiq 2 )] l 5 U a ( Pa ) , (1) 

where a,b,c = 1, 8 and q = p b — p a is the momentum transfer between baryons a and b. 
From Eq. ([I]), the matrix elements of the divergence of the axial currents become 
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b, p b | <PA C | a, p a ) = iU b (p b 



.-(M a + M b )g a A bc ( q 2 ) + q 2 g a 2 bc (q 2 : 



7 5 U a (p a ) , (2) 
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where M a is a baryon mass. Crucial to the derivation of the GTR is the GB pole contribution 
represented in Fig. 1. To explicitly expose the pole term, the matrix element in Eq. (0) can 
be rewritten as 

j\rabc I „2\ 

(b, Pb | 0M< | a, Pa ) = iU b {p b )- f±- l5 U a ( Pa ) , (3) 

\ q z — m A c + it 

where N abc (q 2 ) = g C ab(q 2 )P c (q 2 ) + {q 2 — m 2 )8 abc (q 2 ), m c is a GB mass, and g ca b{q 2 ) the 
baryon-GB form factor, defined such that in the physical basis of the Gell-Mann matrices 
<73,6+«7,6-t7(M^) is equal to g n o nn , etc. P c (q 2 ) represent the couplings of the pseudoscalar 
currents to the GB's, given in the chiral limit by P c = m 2 F c (F c is the decay constant, 
where F n = 92.42 MeV); the q 2 dependence of P c (q 2 ) starts at 0{p i ) and is henceforth 
disregarded. Finally, 5 abc (q 2 ) denotes contributions not involving the GB pole, and it starts 
as a quantity of 0{p 2 ). This separation of pole and non-pole contributions is not unique (the 
off-shell functions separately are not observables) ; for instance, up to higher order terms in 
q 2 , we can choose to remove the q 2 dependence in g ca b(q 2 ) around the point q 2 = m 2 by a 
simple redefinition of 5 abc . 

In the chiral limit d^A ^ = 0, and at q 2 = Eq. (||D gives: 

Mgt(0) = lim q 2 g ab \q 2 ) = F c g cab (0) , (4) 

which is the general form of the GTR. Here M is the common octet baryon mass in the 
chiral limit. In the real world, chiral symmetry is explicitly broken by the quark masses and 
the GB's become massive. In this case, Eqs. (E3) and (Ef) lead to 



> „/„.. >. ,. g ca b(q 2 ) P c (q 2 ) 
m c g 2 ( m c)= } im 2 — yt ~ ■ 5 

q 2 ^ml q z — m/ + le 



In order to define the GTD it is also convenient to take the limit q — > which gives 



(M a + M b )gt(0) = — 2 g C ab{U) P c (0) - 6 abc (0) . (6) 
mi 



The discrepancy A is then defined by 



(M a + M b ) g a A bc (0) = { - ^ gcab ( m 2 c )P c (m 2 c ) . (7) 

m A c 



Notice that while the GTR, Eq. (4), is defined at q 2 = 0, the GTD in Eq. (0) is given 
at q 2 = m 2 because only at that point is the coupling unambiguously determined. At 
leading order in the quark masses, the GTD can then be expressed as follows: 

^ ahc = m 2 c ^- 2 \ogN ahc {q 2 )\ q , =ml . (8) 



II. TREE LEVEL CONTRIBUTIONS 

Throughout we are going to use standard definitions, namely: 

u = exp j , (9) 

X = 2B (s + ip) , (10) 

X± = u ] \u ] ±u\ ] u , (11) 

uj"= ~{u + d ti u-ud>*u + ), (12) 

8* = \l^v v . (13) 

The HBChPT Lagrangian is ordered in powers of momenta and GB masses, which are small 
compared to both the chiral scale and the baryon masses, 

£ = + + £W + . . . . (14) 

Although the Lagrangian is written as a single expansion, it will be useful to keep track of 
the chiral and 1/M suppression factors separately. As will be demonstrated explicitly below, 
leading order (LO) contributions to the GTD appear within £^ 3 \ Subleading contributions 
are suppressed by at least two suppression factors, so we will refer to any contribution at 
the order of £^ as a next-to-leading order (NLO) contribution. 

The tree level contributions to the GTD stem from contact terms in the effective La- 
grangian that can contribute to S abc , and also from terms that can give a q 2 dependence 
to g ca b- First we notice that in HBChPT such terms must contain the spin operator S£ 
that results from the non-relativistic reduction of the baryon pseudoscalar density. There 



are two types of terms which contribute to the GTD. The first type must contain the pseu- 
doscalar source x— The second type must contain monomials such as [D^, u u ]] and 
[D", [Dp, u^]] between the baryon field operators (here T>^ is the chiral covariant derivative). 
However, upon using the classical equations of motion satisfied by the GB fields at 0(p 2 ), 
it turns out that terms of the second type can be recast into terms among which there are 
terms of the first type. In this way, one moves the explicit q 2 dependence from g ca b to contact 
terms, some of which contribute to 5 abc . Such reduction of terms has been implemented for 
C^ 1 ' + C^ 2 ' + £( 3 ', for instance in Ref. |§, and in the relativistic effective Lagrangian as well 
@. Some terms in £( 3 > whose coefficients are determined by reparametrization invariance 
|| may seem at first glance to give a q 2 dependence to g ca b, but a careful calculation shows 
that this is not so. 

Since X- is C(jo 2 ); and since a factor of the spin operator S£ is needed, the LO tree 
contributions to the GTD must come from C^ 3 \ One can further argue that there are no 
contributions from the even-order Lagrangians, C^ 2n \ The reason is that an even number 
of derivatives would require factors in the monomial of the form v .V which, when acting on 
the baryon field, are in effect replaced by V 2 /2M; the other possibility would be factors of 
v.S v that vanish. 

In the case of SU(2), the Lagrangian £( 3 ' has been given by Ecker and Mojzis HI!]. 
There are only two terms in that are of interest to us, namely the terms Oig and O20 
given in Refs. fl8||L0|l. In the scheme used by Ecker and Mojzis these are finite counterterms. 
We note that although Oyj and Oi 8 do contribute to g ca b and to g^ c , they are such that no 



contribution to the GTD results, as noticed in Ref. |10j- In SU(3) there are instead three 
£(3) terms that are of interest to us, namely, 

£<§) rD = iF 19 Tr{BSZ[V,X-, B}) 
+ iD 19 Tt(BSS{V ttX -, B}) 

+ i b 20 Tr(BS»B)Tr(d fM X-)- (15) 
The NLO contributions come from Cqtd an d will not be displayed here. There are, for 



instance, terms quadratic in the quark masses such as Ty(BS£x+[^ /iX-> B]) and others. 
The contribution to 5 abc from C^q TD is given by 

fiabc 2 

ri - 2s [iF l9 f abc + D 19 d abc ] + d cde s d [iF 19 f bea + D 19 d abe ] + s c [-D 19 + b 20 ]6 ab , (16) 



where 

s = ^(m u + m d + m s ) , (17) 
s a = 5 a3 (m u - m d ) - -^=5 a8 (2m s -m u - m d ). (18) 

In deriving Eq. (|16|) from Eq. ( |15|) we used the Ward identity: 

as well as the following correspondence of operators between the heavy baryon and relativis- 
tic theories: 

B v S^pB v ^ -tMB l5P B , (20) 

where B and B v are the relativistic and heavy baryon fields respectively. 

The leading terms in the GTD are therefore of order p 2 . There are several relations among 
the discrepancies that are exact at LO. One of them is the Dashen-Weinstein relation J7|: 



/ \ NNtt 1 / / \ NAK / \ AfSX 

,2 / 9 A \ A NNn 1 2 I o / 9A \ a NAK / 9 A 



A "'"-2<{ 3 {^) A " A ""^J A NZK | . (21 

This particular relation provides useful insight as will be shown in the phenomenological 
discussion. 

Since the bulk of the contribution to the GTD will result from the counterterms of Eq. 
(|T5|), it is important to consider what physics determines their magnitude. It seems likely 
that a meson dominance model may provide the correct picture. In such a model the size 
of the counterterms would be determined by the lightest excited pseudoscalar mesons that 
can attach the pseudoscalar current 575 X a q to the baryons. The relevant such states are 
in the II' octet consisting of 7r(1300), ^(1440) and _fT(1460). The next set of pseudoscalar 



states is in the range of 1800 to 2000 MeV, and thus, one may expect that they only give 
corrections at the order of 20 to 30%. The meson dominance model can be implemented 



using an effective Lagrangian in analogy with Ref. [ Flf . The coupling of the II' octet to the 
pseudoscalar current is obtained from the effective Lagrangian: 

c w = ^(vrrv^n') - 1 -m^ 2 + id^Wx-) + ... , (22) 

where we display only those terms relevant to our problem. Here the IT' octet responds to 
chiral rotations in the same way as the baryon octet. The matrix element of the divergence 
of the axial current is given by 

< I d^Al I W b > = ~d U 'Tr(\ b {\ a ,M q }) 

= -5 ab d Il ,m 2 a , (23) 

and the Il'-baryon coupling can be expressed through the effective Lagrangian: 

C n >B = D'Ti(B l5 {W, B}) + F'Ti(B l5 [W , B}) . (24) 

From Eqs. (E3) and fl2"3| ) one readily obtains the contribution to 5 abc : 

* = -«-A; s «-S' (25) 

Here g$ c B = ^Tr(A b [A c , A a ]) + ^Tr(A 6 {A c , A a }). The current situation is that the couplings 
of the n' are not known, and there is no estimate in the literature that one could judge 
reliable. As we comment later, the GTD's actually serve to determine di[i(q 2 )g^ B much 
more precisely than any model calculation available, provided the meson dominance model 
is realistic. 

III. LOOP CONTRIBUTIONS 

There are several one-loop contributions to the GTD that we illustrate in Fig. 2. There 
are also, at the same NLO, two-loop contributions that we do not display here. Although 
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we do not perform here a full calculation, we do arrive a some interesting observations about 
such NLO effects by loops. Let us consider the loop diagram in Fig. 2a. We can show that 
in HBChPT this loop effect on the GTD is 0(1/M 2 ), and must therefore be suppressed by 
two powers relative to the LO contribution. Indeed, in HBChPT the diagram is proportional 
to the following loop integral: 

w f d d k k^k v 1 + k ■ v/(2M f ) + 0(1/Mf) l + {k + q).v/{2M e ) + 0{l/M 2 ) 



(2n) d k 2 - m 2 k.v + k 2 /{2M f ) - 5m fa (k + q).v + (k + q) 2 /(2M e ) - 6m eb ' 

(26) 



where 5m ab = M a — M^, and T^ u is transverse to the four-velocity v. For spin 1/2 baryons 
in the loop T^ u oc Stfq-S v S". It is also easy to show explicitly that T^ v is transverse if one or 
both lines in the loop are spin 3/2 baryons. From energy-momentum conservation we have 

q , = (M t - M„) - ^. (27) 



Using this and the transversity of T Miy , the expansion of Eq. (|26"D shows no q 2 - dependence at 
(9(1) and 0(1/ M). We conclude that the one- loop diagrams considered here must affect the 
GTD at 0(1/M 2 ) [], and are thus negligible in the large M limit. Another type of one-loop 
contribution is not suppressed by 1/M. These are the diagrams involving the insertion of 
terms from £® as shown in Fig. 2b, which correct the GTD at NLO. Similarly there are 
NLO two- loop contributions that are of leading order in 1/M. 

It is interesting to comment here on a one-loop calculation in the framework of a rela- 
tivistic baryon effective Lagrangian, as used in Refs. P,|T3"||. It turns out that the relativistic 



version of the loop diagram in Fig. 2a gives a finite q 2 dependence to the g ca b coupling, 
namely, 

g cab (q 2 ) - g cab (0) = (-\-) 3 £ g a A fd g e A d 9 f A C J fed (<i 2 , M a , M b , m c ) (28) 



V2F 7 . 

where the integral J^ ed is given by: 



d,e,f=l 



1 For a related discussion, see Ref. IT 
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J fed (q 2 ,M a ,M b ,m c ) 




+ 2(1 - y)M f M b + 2yM b M e + (M f - M e f , (30) 
V(x, y) = (l-x-y) {xM 2 + yM 2 b - M 2 ) - M 2 f x - M 2 y + xyq 2 , (31) 



C(M a ,M b ,M e ,M f ) = (M b + M e )(M a + M f )(M e + M f ) . 



(32) 



One can readily check that for SU(2) one obtains the result in Ref. [0]. 



The interesting thing here is that the contribution to the GTD by the loop is not sup- 
pressed by 1/M. Actually, it is nearly constant for baryon masses ranging from a few 
hundred MeV to an arbitrarily large mass. This result seems at odds with the one from 
HBChPT, but the two can be harmonized as follows: in the limit of large M it turns out 
that in the relativistic calculation there are contributions to the loop integral from momenta 
that are O(M). M acts in fact as a regulator scale. In HBChPT on the other hand, one is 
doing a 1/M expansion of the integrand, which implies that one is assuming a cutoff in the 
loop integrals given by a QCD scale. The relativistic and HBChPT frameworks must each 
lead to the same physical results; in the present case this implies that in order to lead to the 
same results for the discrepancies, the coefficients Fig and -Dig in must be readjusted 
when going from one framework to the other. In the real world, M ~ A x and we may use the 
relativistic calculation as an estimate of this class of loop contributions to the discrepancy. 
For the discrepancies of interest herein, these loop contributions are small, between ten to 
twenty percent of the discrepancies themselves, and smaller than their current errors. The 
numerical results are 



N Nit 
loop 



0.0043 



(33) 



NAK 
loop 



-0.044 



(34) 



NT,K 
loop 



0.044, 



(35) 



where we use D = 0.79 and F = 0.46 for the SU(3) axial vector couplings. 



8 



Of course the calculated loop contribution is not all that there is; the inclusion of decuplet 
baryons in the loop also gives contributions to the discrepancy. (Ref. |TJj] discusses some 
A(1232) effects with only two quark flavors.) Using Rarita-Schwinger propagators and three 
quark flavors, we have checked that the g 2 -dependent part does show an UV divergence in 
the relativistic framework. HBChPT also permits two-loop contributions at NLO. Currently 
a more complete calculation of the discrepancy at NLO is underway. [15] 



IV. RESULTS 

There are only three discrepancies that can be determined from existing data on baryon- 
pseudoscalar couplings: A NNn , A NAK , and A NT,K . 

Due to the smallness of the u and d quark masses, A NNn is necessarily very small, and 
its determination requires a very precise knowledge of the g^NN coupling (g^ and F n are 
already known to enough precision, leaving most of the uncertainty in the determination 
of A NNn to the uncertainty in g^NN)- The most recent determination of g n NN from NN, 
NN and irN data is by the Nijmegen group |TjJ. They analyzed a total of twelve thousand 
data and arrived at g n ^N = 13.05 ± 0.08. Similar results are obtained by the VPI group 
|17| . Since the errors quoted are only statistical, in our fit below we will increase the error 
by about a factor of two in order to roughly account for systematic uncertainties. There is 
still some disagreement between determinations of g n NN by different groups. Larger values 
have been obtained, such as g^NN — 13.65 ± 0.30 by Bugg and Machleidt fL8| , and a similar 



result by Loiseau et al. ||19|| . As we find out below, our analysis of the discrepancies strongly 
favors the smaller g wNN values. Using F w = 92.42 MeV, i^j = -1.267 ± 0.004 0, 
Eq. (0) gives, 

A S W = °- 014 ± °- 006 for 9*nn = 13.05 ± 0.08, (36) 
A S W = 0.056 ± 0.020 for g wNN = 13.65 ± 0.30. (37) 
The determination of the gxNA and gxNT, couplings relies on a more sparse data set. 



The Nijmegen group analyzed data from YY production at LEAR, and they obtained |2T 



9kna = 13.7 ±0.4 and gxNT. = 3.9 ±0.7. These values are consistent with an earlier analysis 



by Martin |EJ, where only an upper bound for Qknt, is given. Using = 1.22 F n and 



/ \ NAK / \ NT,K , „ _ 

(Mj = -0.718 ± 0.015 and (^-J = 0.340 ± 0.017 |2(|, Eq. (0) gives, 



Af x Jf = 0.17 ±0.03 (38) 



A™ = 0.17 ±0.14 (39) 

Disregarding SU(2) breaking, which implies that there is no contribution from the term 
proportional to 6 20 to these discrepancies, we can use the three measured discrepancies to 
determine the two LO parameters in HBChPT: 

MF 19 = 0.4 ±0.1 GeV" 1 , (40) 
M Di 9 = 0.7 ± 0.2 GeV -1 , (41) 

where M is here the common baryon-octet mass in the chiral limit. Both choices for g^NN, 
given in Eqs. ([!6|) and (0), lead to values for M Fig and M Dig that agree within the 
quoted uncertainties. The LO discrepancies resulting from our fit are: 

A NNn = 0.017; 0.018, (42) 
A™ = 0.17; 0.18, (43) 
A NJ]K = 0.17; 0.19, (44) 

where the quoted results correspond respectively to the smaller and larger g^NN couplings. 
The larger value A^"" = 0.056 of Eq. (|37D, corresponding to the larger g n NN coupling, 
cannot come out consistently from the fit. To understand this one can use the Dashen- 
Weinstein relation, Eq. (pi]), which holds exactly in our LO calculation. For the results 
of the discrepancies involving the hyperons the term proportional to A NT,K in the Dashen- 
Weinstein relation is about one fifth of that proportional to A NAK , and the right hand side 
of Eq. (|2T| ) would imply that A NNw must be about 1.5%. The only way to accomodate a 
larger A NN?T would be larger A NAK and A NT;K or else a large deviation from the Dashen- 
Weinstein relation. The latter seems unlikely because the corrections to the relation must be 

10 



suppressed by two powers in HBChPT (this is so because the corrections to the axial- vector 
couplings and to the discrepancies are of 0{p 2 )). On the other hand the former possibility 
would require that the magnitudes of Qknk and gxNT, be unrealistically large. In fact, A NAK 
and A NT,K should be close to unity, implying a serious failure of the low energy expansion. 
Thus, we conclude that only the smaller values of A NNw , and thus of g w NN, are consistent. 
This shows the importance of the current analysis of the GTD in SU(3). 

Finally, the coupling constants required in the meson dominance model resulting from 
our analysis are as follows: 

dn>F' = 2.4 ± 0.5GeV (45) 
d n >D' = 4.5 ± 0.5GeV. (46) 

Since here F' and D' are baryon-meson couplings, it is not unreasonable that they should 
have values similar to those of, say, the pion-nucleon coupling. This would imply that the 
coupling dw should be a few hundred MeV. This makes the meson dominance picture quite 
plausible. 

In conclusion, we have shown that the GTD in SU(3) is given at leading order by two 
tree-level contributions, and that the corrections are suppressed by two powers in HBChPT. 
Some of the loop corrections were calculated explicitly and found to be small. Our leading 
order analysis indicates a strong preference for a smaller Goldberger-Treiman discrepancy 
in the pion-nucleon sector, thus favoring the smaller values of the pion-nucleon coupling 
extracted in recent partial wave analyses. 
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FIGURES 




FIG. 1. Diagrams representing the contact term and the pole term in the matrix elements of 
the divergence of the axial currents. Crosses represent the divergence of axial currents. 
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FIG. 2. One loop diagrams that give NLO corrections to the GTD. In (a) the cross represents 
the divergence of the axial current obtained from the 0(p 2 ) Lagrangian, and in (b) the same 
divergence obtained from the 0(p s ) baryon Lagrangian. 
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